Quantum Phase Transition in the Sub-Ohmic Spin-Boson Model: Extended 

Coherent-state Approach 
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We propose a general extended coherent state approach to the qubit (or fermion) and multi- 
mode boson coupling systems. The application to the spin-boson model with the discretization of 
a bosonic bath with arbitrary continuous spectral density is described in detail, and very accurate 
solutions can be obtained. The quantum phase transition in the nontrivial sub-Ohmic case can be 
located by the fidelity and the order-parameter critical exponents for the bath exponents s < 1/2 
can be correctly given by the fidelity susceptibility, demonstrating the strength of the approach. 



PACS numbers: 05.10.-a, 03.65.Yz, 71.10.Fd, 05.30.Jp 

The spin-boson model [J describes a qubit (two- 
level system) interacting with an infinite collection of 
harmonic oscillators that models the environment act- 
ing as a dissipative bosonic " bath" . There are currently 
considerable interests in this quantum many-body system 
due to the rich physics of quantum criticality and deco- 
herence , applied to the emerging field of quantum 
computations, quantum devices @, and biology @, Q. 
The dissipative environment Q in the spin-boson model 
is characterized by the spectral function J(u>) with fre- 
quency behavior J{u>) oc lo s . The spin-boson model un- 
dergoes a second-order quantum phase transitions (QPT) 
from delocalized to localized phase with a sub-Ohmic 
bath (0 < s < 1) and a Kosterlitz-Thouless type transi- 
tion in the Ohmic case (s = 1). 

To provide reliable solutions for the spin-boson model, 
a typical multi-mode system, is quite challenging. Al- 
though several numerical methods applied to the sub- 
Ohmic case EH, EH can reproduce the phase 
diagram, only recent quantum Monte Carlo (QMC) sim- 
ulations [l5| and exact-diagonalization studies [l6| are 
capable of correctly extracting the critical exponents in 
the QPT. The critical behavior in the previous stan- 
dard numerical renormalization group (NRG) calcula- 
tions [l(| EH, EII 1S incompatible with a mean-field tran- 
sition due to the failure of the quantum-to-classical map- 
ping with long-range interactions for s < i. More re- 
cently, the early standard NRG results were improved by 
a modified NRG algorithm [l7j. and the mean- field be- 
havior for s < 1/2 was also reproduced. 

In this paper, we present a general accurate approach 
to the qubit (or fermion) and multi-mode boson coupling 
systems. As a important example, we focus on the sub- 
ohmic spin-boson model here. It can also be easily ex- 
tended to the famous Holstein model [HI and the multi- 
mode Dicke model [20j. The crucial procedure is to em- 
ploy extended coherent states to represent the bosonic 
states. The QPT in the sub-ohmic spin-boson model will 
be analyzed by means of the quantum information tools, 



such as the ground state fidelity and fidelity susceptibil- 
ity [2TM23I ] . It is a great advantage to use the fidelity to 
characterize the QPT, since there should be a dramatic 
change in the fidelity across the critical points. More- 
over, the non-trivial order-parameter critical exponents 
can be obtained with scaling of the fidelity susceptibility. 
The Hamiltonian of the spin-boson model is given by 
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where a x and a z are Pauli matrices, A is the tunneling 
amplitude between two levels, u> n and a) n are the fre- 
quency and creation operator of the n-th harmonic oscil- 
lator, and A n is the coupling strength between the n-th 
oscillator and the local spin. The spin-boson coupling is 
characterized by the spectral function, 

J{lo) = A 2 t (5(w„ -u)) = 2irauj 1 c ~ s uj s , < us < uj c 

n 

(2) 

with w c a cutoff frequency. The dimensionless parameter 
a denotes the strength of the dissipation, s = 1 stands 
for an Ohmic dissipation bath. The rich physics of the 
quantum dissipation is second-order QPT from dereal- 
ization to localization for < s < 1, as a consequence of 
the competition between the amplitude of tunneling of 
the spin and the effect of the dissipative bath. 

We here propose a solution of the spin-boson model 
by exact diagonalization in the coherent-states space. To 
implement our approach, we first perform discretization 
of the bath speciation function, according to the loga- 
rithmic discretization of the continuous spectral density 
J {(J) in the NRG [13, EH, El • The discrete Hamiltonian 
is therefore expressed as 



n * n 

(3) 
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with 

6i=7,7 2 / dxJ(x)x,j% = / da;J(x) 

In order to ensure the convergence of the results, the 
disctetiztion parameter is chosen A = 2. 

The present basic scheme is similar to that in the 
single-mode Dicke model [25| and the two-site Holstein- 
Hubbard model [26j]. For convenience, we assume that 
| ipi) and 1^2) are the bosonic states corresponding to 
spin up and down. Introducing a displacement shift pa- 
rameter g„ — 2 £ l 'l/¥ 13' we propose the following two 
coherent bosonic operators 



A„ — a n + 9n i A n ~ a n + <7ri 

-6,1 = a£ - g n .B n = a n — g„ 



(5) 
(6) 



The corresponding vacuum states |0)^ n and |0)b„ are 
just the coherent states in a n with eigenvalues ~^g n in 
terms of e T9 ™° ~ 9 "^ 2 |0) ajl . \n k )A k and \n k )B k correspond 
to Fock states of the new bosonic operators and B^ 
with n k bosons for a frequency £Jfe. and \ip2) can be 
expanded in the bosonic coherent states of a series of rife, 
which arc orthonormalized in the new bosonic operators 
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where c',„ x are coefficients with respect to a series of 
{ni 1 n2-..n]\[} for different bosonic modes, and N tr is the 
bosonic truncated number. 

Then the Schrodinger equations of the Hamiltonian (3) 
are derived as 
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After the substitution of Eqs. (7) and (8) and Left mul- 
tiplying the bosonic coherent states with the both sides 
of Eqs. and ([JTJ, we have 
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The bosons state |rifc) and \rrik) with different coherent 
bosonic operators and B^ are not orthogonal. The 
overlap can be denoted by A k {m k \n k ) Bk = {-l) nk D mkTlk 
and B k {m k \n k ) Ak = {~l) mk D mkTlk with 
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According to the symmetry of Hamiltonian for e = 0, 
the coefficients satisfy c mi ... mN = ±(-l)^<= " lfc cZ mi ... mjv . 
Eqs.( ITT]) and ( [T2")) can then be transformed into the 
following set of coupled equations 

A 00 

2J c i^} n ^"fc+E^^"^ 2 ) !™^ = Sc {™fc} 

«i...njv fe=l i=0 

(13) 

A complete implementation of the numerical diago- 
nalization is described below to obtain the amplitudes 
set of cs nk \ of the bosonic state ipi (if2 )■ The Hilbert 
space can be labeled by a vector n = (n%, ...,npf) with 
n k — 0, 1, ...,Ntr- The sum of bosonic number n k is re- 
stricted to truncated number Ntr, e.g. J2 n k < N tr - For 
example, with a set of A~ — 3, Aft r = 3 the involved con- 
figurations of bosonic states ri2, n>3 > are expressed 
as following: 

|000), 
|100), |010), |001), 
|200),|110), |020), |101), |011), |002), 

|300), |210), |120), |030), |201), |111), |021), |102), |012), |003). 

Consequently, the total number of basis states N s = 20. 
To obtain the true exact results, in principle, the num- 
ber of bosonic modes A" and the truncated number N tr 
should be taken to infinity. Fortunately, in the present 
calculation, setting A~ = 16, which is big enough in NRG, 
and N tr — 5 is sufficient to give very accurate results 
with relative errors less than 10~ 5 in the whole param- 
eter space. The following results are just obtained with 
N = 16 and N tr = 5. 

To study the QPT in the sub-Ohmic spin-boson model, 
we employ the ground state fidelity to locate the critical 
point a c . A simple expression of the ground-state fidelity 
is given just by the modulus of the overlap 

F(a,a') = \{ Vl (a)\ Vl (a')) + (^ 2 (a)|^ 2 (a'))| (14) 

The QPT is expected to be signaled by a drop in 
the fidelity corresponding to two arbitrarily neighboring 
Hamiltonian parameters a' = a + 5a (2~ll. l23j] . Based on 
the normalized ground states \<pi) and \f2), we now illus- 
trate our results obtained by numerically diagonalization 
ofEq.([ni). 
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Fig. Q] (a) shows the behavior the fidelity F(a,a') 
in the sub-Ohmic case with s = 0.5, 0.6, 0.8, 0.9 for the 
spin tunneling amplitude A = 0.01 . A sharp drop at 
the critical point a c separates the delocalized phase at 
small a and the localized phase at large a. So it is evi- 
dent that we can locate the critical points a c efficiently 
by the ground state fidelity. It is interesting that the 
ground state fidelity does not drop to at critical points, 
demonstrating a continuous QPT[24|. 

The QPT from delocalized to localized phases can also 
be shown by behavior of the tunneling motion (a x ) be- 
tween spin up and spin down (l3j |. (<r x ) goes rapidly to 
zero in the localized phase and is finite in the delocal- 
ized phase. The a(A, s) dependence of (a x ) is shown 
in Fig. [IJb). For all value of s, we observe that (a x ) is 
continuous at the transition. The discontinuous behavior 
observed previously [l3[ may be attributed to the special 
variational approach itself, and is perhaps worthy of a 
further study. 

As discussed above, both the fidelity F(a, a') and tun- 
neling parameter (a x ) can be used to locate the critical 
points of the QPT. We observe that both quantities can 
give nearly the same critical points. The phase bound- 
aries obtained by either methods as a function of s is 
plotted in Fig. [5] in case of the tunnel splitting A ranging 
from 10~ 4 to 10 _1 for e = 0. The results from previous 
NRG techniques [T3, E| are also collected for compari- 
son. It is interesting to note that the present results for 
the critical points are in good agreement with the NRG 
ones. Because we also use the truncated NRG Hamilto- 
nian (3), the critical points should be slightly above the 
QMC ones[l5| where all frequencies are included. For 
fixed truncated Hamiltonian , as N tr increases, a c con- 
verges to the true value from above very quickly. We 
believe that we obtain the converging critical points for 
any values of A for fixed N in the present work. 

The ground-state fidelity susceptibility y is defined as 
the second derivative of the fidelity [U |23| | 
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FIG. 1: (Color online), (a) The fidelity F as a function of a 
for various values of s and (b) The tunneling (a x ) between two 
states of the spin as a function of a in the case of A = 0.01 
and e = 0.. 




x(a) = 2 lim 

5a— *0 



1 - F(a, 8a) 
So 2 



(15) 



Since x is independent of the arbitrary small parameter 
Sa , it is regarded as a more effective tool to detect the 
singularity in QPT. As addressed in Ref. [231 ]. the fidelity 
susceptibility is similar to the magnetic susceptibility. In 
the localized phase, the scaling behavior of fidelity sus- 
ceptibility x obeys 



X(a) oc \a - 



(16) 



Recently, the studies from the QMC approach [15j, the 
exact-diagonalization [16| . and the modified NRG [TtJ 
have shown that the quantum-to-classical mapping is 
valid for the sub-Ohmic spin-boson model, i.e. the 
critical exponents are classical, mean-field like, in con- 
trast with the early standard NRG calculations where 
the quantum-to-classical mapping is suggested to fail for 
s < 1/2. It was argued [l5| that the standard NRG is 



FIG. 2: (Color online). The delocalized-localized transition 
point olc as functions of s obtained by the present approach 
in the case A = 10~\ 10~ 2 , 10" 3 , 10" 4 . The NRG data are 
also shown for comparison. 



not able to capture the correct physics in the localized 
phase. It is known that the localized phase is two-fold 
degenerate. Our ansatz Eqs. (7-10) is just proposed for 
these two states, and the unknown coefficients can be 
obtained by solving the Schrodinger equations very accu- 
rately. Therefor, we will extract the susceptibility critical 
exponent for s < 1/2, to address this crucial controversy. 

Fig. [3] presents that the fidelity susceptibility x as a 
function (a — a c ) for s = 0.1, 0.2, 0.3, 0.4 in log-log scale. 
All curves show almost perfect straight line with a slope 
very close to —1, demonstrating that the susceptibility 
critical exponent may be just equals to (3 = —1 with 
X ~ (a — a c )~ l . Recently, the critical exponent of the 
magnetic susceptibility has been estimated to be —1 by 
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FIG. 3: (Color online). The scaling behavior of the fi- 
delity susceptibility x ~ \ a ~ a c\ 13 in log-log scale for s — 
0.1,0.2,0.3,0.4 with A = 0.1. The data fit well with the 
straight line with slope very close to —1. 

QMC simulations 15] and exact-diagonalization stud- 
ies Hal. We do not think this is a coincidence. You 
et al [23[ have shown a neat connection between the fi- 
delity susceptibility \ an d the magnetic susceptibility \m 
through x = Xm/^KsT (Kb is the Boltzmann constant, 
T is the temperature). We believe this relation is also 
applicable to zero temperature, and therefore these two 
susceptibilities can give the same order-parameter crit- 



ical exponents in QPT. We also confirm the mean-field 
behavior for 1/2 < s < 1. 

In summary we have introduced an efficient algorithm 
in the new bosonic coherent Hilbert space and presented 
reliable solution for the sub-Ohmic spin-boson model. 
The ground-state fidelity, which is a quantum informa- 
tion tool, is employed to locate the critical coupling 
strength a c of the QPT . The transition from the lo- 
calized phase to delocalized phase is accompanied by a 
minimum of the fidelity. Furthermore, the fidelity sus- 
ceptibility gives the order-parameter critical exponent 
P = — 1 in the case s < 1/2, which agrees well with 
the exponent of magnetic susceptibility. Both behaviors 
of the tunneling (a x ) and the fidelity around the criti- 
cal point exclude the possibility of the first-order QPT. 
We stress that all eigenstates and eigenvalues of the spin- 
boson model can be obtained accurately and many ob- 
servable can be calculated directly within the present ap- 
proach. The present technique to deal with bosons would 
be combined with other established methods. 
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